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The iron pnictide and chalcogenide compounds 
are a subject of intensive investigations due to 
their high temperature superconductivity. [Ij They 
all share the same structure, but there is significant 
variation in their physical properties, such as mag- 
netic ordered moments, effective masses, supercon- 
ducting gaps and Tc. Many theoretical techniques 
have been applied to individual compounds but no 
consistent description of the trends is available [2j. 
We carry out a comparative theoretical study of a 
large number of iron-based compounds in both their 
magnetic and paramagnetic states. We show that 
the nature of both states is well described by our 
method and the trends in all the calculated physical 
properties such as the ordered moments, effective 
masses and Fermi surfaces are in good agreement 
with experiments across the compounds. The varia- 
tion of these properties can be traced to variations in 
the key structural parameters, rather than changes 
in the screening of the Coulomb interactions. Our 
results provide a natural explanation of the strongly 
Fermi surface dependent superconducting gaps ob- 
served in experiments [3] . We propose a specific op- 
timization of the crystal structure to look for higher 
Tc superconductors. 

The iron pnictides are Hund's metals [4 , where the inter- 
action between the electrons is not strong enough to fully 
localize them, but it significantly slows them down, so that 
the low energy quasiparticles have much enhanced mass. 
These quasiparticles are composites of charge and a fiuc- 
tuating magnetic moment originating in the Hund's rule 
interactions which tend to align electrons with the same 
spin and different orbital quantum numbers when they find 
themselves on the same iron atom. 

A central puzzle in this field is posed by the varia- 
tion of the ordered magnetic moment across the iron pnic- 
tides/chalcogenides series. In the fully localized picture the 
atom resides in a single valence, therefore the ordered mo- 
ment is equal to the atomic moment per iron), pos- 
sibly reduced by quantum fiuctuations. This picture is re- 
alized in cuprate superconductors where quantum fiuctua- 
tions reduce the Cu^+ moment by 20%. In the fully itinerant 
weak coupling picture, such as spin density wave (SDW) in 
chromium metal, the ordered moment is related to the de- 
gree of Fermi surface nesting. It is by now clear that the iron 
pnictides are not well described by either fully localized or 
fully itinerant picture, nor by the density functional theory 
(DFT), which greatly overestimates the ordered magnetic 
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FIG. 1: Ordered magnetic moments and mass enhance- 
ments in iron-based compounds, (a)the DFT+DMFT cal- 
culated and experimental 6 - 13 Fe magnetic moments in the 
SDW and DSDW states. Also shown is the calculated fluctuating 
moment in the PM state. (b)The DFT+DMFT calculated mass 
enhancement m*/mband of the Fe 3d orbitals in the PM state and 
the low energy effective mass enhancement obtained from opti- 
cal spectroscopy experiments p'MTS] and (angle-resolved) pho- 
toemission spectroscopy experiments p!9H24] . 



moments. It has been advocated that the shortcomings 
of DFT can be circumvented by incorporating the physics 
of long wavelength ffuctuations [5]. Here we take the op- 
posite perspective. While critical long- wavelength ffuctua- 
tions certainly play a role near the phase transition lines, 
we will show that the local ffuctuations on the iron atom 
can account for the correct trend of magnetic moments and 
correlation strength in iron pnictides/chalcogenide layered 
compounds. 

Using the combination of density functional theory and 
dynamical mean ffeld theory (DFT+DMFT) (see online ma- 
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terial for details), we studied two different real space or- 
derings, the SDW ordering, characterized by wave vector 
(7r,0,7r) (this vector is written in coordinates with one Fe 
atom per unit cell), which is experimentally found in iron ar- 
senide compounds, and (7r/2, 7r/2, tt) ordering, denoted by 
double stripe spin density wave (DSDW). The latter was 
found experimentally in FeTe. Figure [TJa) shows our the- 
oretical results for the ordered moment in both phases to- 
gether with experimentally determined values [6HT3] from 
across all known families of iron-based superconducting 
compounds. There is an overall good agreement between 
theory and experiment, in particular LaFePO is predicted 
to be nonmagnetic, the majority of 1111 and 122 compounds 
have ordered moment below 1.0 /iB, and FeTe orders with 
DSDW moment of 2.1 /i^- 

We now explain the variation of the ordered moment in 
terms of real space and momentum space concepts. The 
size of the fluctuating local moment, which can be extracted 
from neutron scattering experiments, gives an upper bound 
to the size of the ordered magnetic moment and is also plot- 
ted in Fig. [ija). It was computed from the atomic histogram 
displayed in Figj2|c), which shows the percentage of time 
the iron 3d electrons spend in various atomic configurations 
when the system is still in its paramagnetic (PM) state. 
Only high spin states, which carry a large weight as a result 
of the Hund's rule coupling in iron, are displayed (see also 
online material for complete histogram). A monovalent his- 
togram with only the atomic ground state would give iron 
magnetic moment of 4 /i^. 

In a correlated Fermi liquid, the spin excitations are de- 
scribed in terms of individual particle hole pair excitations 
and their collective motion. Their residual interaction can 
lead to a magnetic state when the particle hole excitations 
condense at non-zero wave vector. A large quasiparticle 
mass, naturally facilitates this condensation, hence we ex- 
pect that the size of the ordered moment will correlate with 
the mass of the quasiparticles. In figure [TJb) we display 
separately the quasiparticle mass for all iron 3d orbitals 
in the PM state and we normalize it to its band value. 
Clearly there is some correlation between mass enhance- 
ment in Figjljb) and ordered magnetic moment in Figjlja) 
across various families of iron-based compounds. In par- 
ticular, correlations are too weak for ordering in LaFePO, 
while very heavy quasiparticles in FeTe produce large mo- 
ment of 2.1 jiB' However, there are other factors presented 
below, such as kinetic frustration, orbital differentiation and 
Fermi surface shape, which together conspire to produce the 
magnetic orderings displayed in Figjlja). 

The quasiparticle mass displayed in Fig. [ijb) is quite 
moderate in phosphorus 1111 compound on the right hand 
side of Figjljb), but correlations are significantly enhanced 
in arsenic 122 and 1111 compounds. Notice however, that 
enhancement is not equal in all orbitals, but it is signif- 
icantly stronger in the t2g orbitals, i.e., xz^ yz^ and xy. 
The correlations get even stronger in 111 compounds, such 
as LiFeAs and NaFeAs, and finally jump to significantly 
larger values of the order of five in selenides KFe2Se2 and 
CsFe2Se2. Finally, the mass enhancement of the xy orbital 



in FeTe exceeds factor of seven compared to the band mass, 
which is typical for heavy fermion materials, but is rarely 
found in transition metal compounds. We displayed only 
a lower bound for this mass as the end point of an ar- 
row in Fig. [TJb), because the quasiparticles are not yet well 
formed at studied temperature T=116K. Notice the strong 
orbital differentiation in FeTe, with xz/yz mass of five and 
eg mass enhancement of only three. This orbital differenti- 
ation signals that the material is in the vicinity of an orbital 
selective Mott transition, as proposed previously for other 
iron pnictides [14^, where xy orbital is effectively insulat- 
ing while other orbitals remain metallic. In Figjllb) we also 
display mass enhancement extracted from optics (TSHIH] and 
ARPES[19 - 24 measurements, and notice a good agreement 
between our theory and experiment when available. The ef- 
fective mass extracted from ARPES and optics should be 
compared with that of the t2g orbitals which contribute 
most of the spectral weight at low energy. 




FIG. 2: Structure, orbital occupation and probability of 
selected atomic states of Fe. (a)The Fe-X (X—P, As, Se and 
Te) distance and X-Fe-X angle in iron-based compounds, where 
the two X atoms are in the same ab plane. Note this angle is 
different from the X-Fe-X angle where the two X atoms are in 
different ab planes; (b)The orbital occupation of the xy orbital 
and the average values for the eg orbitals and all five orbitals; 
(c)The probability of selected atomic configurations of Fe where 
N (S) is the total number (spin) of Fe 3d electrons in the atomic 
configuration. 

The large mass enhancement in Hund's metals is due to 
an orbital blocking mechanism. If the Hund's coupling is 
very large, only the high spin states have a finite probabil- 
ity in atomic histogram. The atomic high spin ground state 
has maximum possible spin 5* = 2, and is orbit ally a sin- 
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glet, which does not ahow mixing of the orbit als and leads to 
orbital blocking, i.e., {gs\dl^df3\gs) = when a ^ /3, where 
1^5) is the atomic ground state in the 3d^ configuration, and 
a is the iron orbital index. In the localized limit and in the 
absence of crystal field effects, it is possible to derive a low 
energy effective Kondo model, which has Kondo coupling 
for factor of (2S' + 1)^ smaller than a model without Hund's 
coupling [25] . Since the Kondo temperature Tk depends on 
the Kondo coupling /q exponentially {Tk cx exp(— I/Jq)), 
this results in enormous mass enhancement of the order of 
exp(((25' + 1)^ — l)//o) compared to the system with negli- 
gible Hund's coupling (see also online material). 

Having established why heavy quasiparticles form in iron 
pnictides and chalcogenides, we can now study how the key 
parameters of the crystal structure control the strength of 
correlations and other physical properties, keeping the same 
on-site Coulomb interaction matrix. The iron-pnictogen dis- 
tance, displayed in Fig. |2|a), controls the overlap between 
iron and pnictogen atom and hence makes iron electrons 
more localized (itinerant) with increasing (decreasing) dis- 
tance. The largest distance is achieved in compounds with 
larger chalcogenide ion, such as in FeTe, which results in 
very heavy quasiparticles, as seen in Fig. [ijb). The vari- 
ation in distance alters the overall bandwidth moderately. 
The Hund's orbital blocking mechanism amplifies this vari- 
ation. 

The second key structural parameter is the tetrahedron 
shape, which is parameterized in terms of pnictogen-Fe- 
pnictogen angle, displayed in Fig. [2|a). This angle is equal 
to 109.5° for an ideal tetrahedron, and is much smaller in 
FeTe, where the Te ion is pushed further away from Fe 
plane. The shape of the tetrahedron controls the crystal 
field levels, which in turn control the orbital occupancies. 
We display them in Fig. [2|b). The average occupation of 
iron atom is around Ud = 6.25 across all the compounds 
studied, which leads to an average orbital occupation of 
Ua = 1.25. A deviation from ideal angle enhances the crys- 
tal field splittings between xy and the degenerate xz/yz 
orbital and also changes the splitting between t2g and eg 
orbitals. Heavier quasiparticles with smaller quasiparticle 
bandwidth are more susceptible to the crystal field split- 
ting, hence the orbital differentiation is largest in FeTe but 
very small in LaFePO. The net result of crystal field split- 
tings and quasiparticle mass is the charge transfer from the 
t2g to eg orbitals as seen in Fig. [2|b), and among t2^'s 
the xy orbital loses most charge with increased correlation 
strength, pushing its occupancy closer to integer filling. 

Furthermore, the effective hopping between neighboring 
iron atoms has two contributions, one is due to direct iron- 
iron overlap, and the second is indirect hop through pnicto- 
gen atom. The two contributions to the diagonal hopping 
ta,a have opposite sign and destructively interfere. For the 
xz and yz orbital, the indirect hopping through pnictogen 
is larger than direct iron-iron hop. For the xy orbital, the 
two contributions are very similar, and when the pnictogen 
height is sufficiently large, such as in FeTe, the indirect hop 
is reduced and the two contributions almost exactly cancel 
each other, resulting in vanishing effective nearest neighbor 



iron-iron txy^xy hopping. This kinetic frustration mechanism 
contributes to the dramatic enhancement of the xy mass in 
the FeTe compound. 
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FIG. 3: Fermi surface. The DFT+DMFT (top row) and 
DFT (bottom row) calculated 2D Fermi surface in the V plane 
[kz = 0) for (a)LaFePO; (b)BaFe2As2; (c)LiFeAs; (d)KFe2As2. 
The Fermi surface is colored in red, green and blue according to 
its orbital character of xy, xz and yz^ respectively. 

In itinerant systems, the shape of the Fermi surface, or 
the Fermi surface nesting is relevant for deciding which mag- 
netic ordering wave vector is realized when the residual in- 
teractions among the quasiparticles is sufficiently strong. 
Additional terms, arising from the incoherent part of the 
electron become increasingly important as the localization 
threshold is approached. 

In Fig. [3] we display DFT+DMFT Fermi surface together 
with DFT predictions. In the moderately correlated end, 
such as in the phosphorus 1111 compounds, our theoretical 
predictions match DFT results. However, when correlations 
become sizable, such as in LaFeAsO or BaFe2As2, the xy 
orbital starts to play a special role, which results in slightly 
modified Fermi surface shape and character compared with 
DFT, while respecting the Luttinger theorem. In BaFe2As2 
DFT predicts that the outer pocket at F is of xz/yz char- 
acter, while DFT+DMFT predicts that the outer pocket is 
of xy character, in agreement with experiments [26]. This 
effect of growing xy pocket at F and consequently shrinking 
of xz/yz pocket is even more apparent in LiFeAs. In the lat- 
ter compound DFT also predicts the outer pocket to be of 
xy character, but its size is considerably smaller than mea- 
sured in ARPES experiment [21 . DFT+DMFT increased 
size of the outer xy pocket together with the butterffy shape 
of the xz/yz pocket at F are in good agreement with exper- 
iment of Ref.21, These changes in the shape of the Fermi 
surfaces are the momentum space counterpart of the real 
space picture of charge transfer among the iron Zd orbitals 
shown in Fig. [2|b). This is because the decrease (increase) 
of the xy {xz^yz) orbital occupancy results in the increase 
(decrease) of the hole pocket size. Finally, the Fermi sur- 
face of KFe2 As2 displayed in Fig. [S] has only hole pockets 
around F but no electron pockets at M, hence there is no 
Fermi surface nesting to facilitate the long range magnetic 
order. Indeed KFe2As2 can not sustain SDW ordering and 
only a tiny DSDW moment can be stabilized, as shown in 
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Fig. [TJa). Even though the mass enhancement in KFe2As2 
and 111 compounds is substantial, the Fermi surface nesting 
stih plays an important role in stabilizing magnetic order- 
ing. 

The fluctuating moment presented in Fig.jlja) monotoni- 
cally increases with increased correlation strength, and con- 
stitutes an upper bound to the size of the ordered magnetic 
moment. However, even when the Fermi surface nesting is 
quite good, such as in 1111 and many 122 compounds, the 
ordered moment is substantially reduced from this upper 
bound. A part of the reduction is due to kinetic frustra- 
tion, discussed above. This effect is properly treated by 
DFT method, nevertheless the ordered moment predicted 
by DFT is around 2.0 /i^ (see also online material), sub- 
stantially above the experimentally measured values. In 
the DFT+DMFT theoretical method, the orbital differenti- 
ation is responsible for large overall reduction of the static 
moment. In very itinerant system, such as LaFePO, the 
quasiparticles are too weakly interacting to condense, hence 
moderate correlations with mass enhancement of 1.5 do not 
sufficiently localize electrons to allow magnetic ordering. In 
most of other compounds, the localization and hence the 
effective mass increase is substantial only in t2g orbit als, 
while eg orbitals remain only moderately correlated. Conse- 
quently the ordered magnetic moment is small in orbitals, 
which causes only a fraction of the fluctuating moment to 
order. Only when correlations are very strong and equal 
in all orbitals, almost the entire fluctuating moment orders. 
Such an example is provided by K0.8Fe1.6Se2, which is ob- 
tained by introducing iron vacancies into KFe2Se2, where 
the entire fluctuating moment of 3.3/iB orders. |23 

We conclude with some experimental consequences of the 
theory. We established that in compounds with substantial 
mass enhancement, the xy orbital is the heaviest and most 
incoherent, placing FeTe at the verge of an orbitally selective 
Mott transition. Furthermore we have identified the chem- 
ical handle, kinetic frustration, responsible for this effect. 
This idea can be tested by applying uniaxial pressure on 
the FeTe, which should result in a noticeable restoration of 
coherence in the transport properties. Our results also sug- 
gest a natural origin for a particle-hole asymmetry in doping 
the parent compounds. Reducing the iron occupancy of the 
3d orbital, brings the occupancy of the xy orbital closer to 
unity, and increases the correlation strength, which in turn 
strengthens the magnetic moment. This has been observed 
in ARPES studies of the BaFe2As2 family [28j. 

In the magnetic state, the in-plane resistivity is very 
anisotropic, as has been shown in optical and transport 
studies of BaFe2As2 [29l [30]. This is a consequence of 
the strong low energy orbital polarization of the xz and 
yz orbital [30 . Increased correlation strength results in an 
increasing participation of the xy orbital, with magnetic mo- 
ment concentrated in the xy orbital, which does not cause 
the in-plane polarization. Hence, in-plane transport and 
optical anisotropy should be reduced with increase of corre- 
lation strength and orbital differentiation. 

Our work suggests that larger mass enhancement of the 
xy orbital leads to smaller superconducting gap on the most 



outside hole pocket centered at F, which is mostly of xy 
character. Furthermore, large degeneracy is a fertile ground 
for superconductivity, while large orbital differentiation is 
harmful, which suggests that superconductivity is hard to 
achieve when X-Fe-X angle is small. On the other hand, 
large X-Fe-X angle in iron pnictides and chalcogenides is 
accompanied by small X-Fe distance, which weakens mag- 
netism and hence likely undermines superconducting pair- 
ing strength. We thus suggest that good candidates for high 
temperature superconductivity are compounds with X-Fe- 
X angle close to ideal angle, as observed by Lee et a/., [31, in 
order to achieve small orbital differentiation, but with the 
largest possible X-Fe distance to strengthen spin fluctua- 
tions. 

METHOD 

We use fixed Coulomb interaction parameters for all ma- 
terials in our DFT+DMFT calculations in order to keep 
parameter-free spirit and to demonstrate that the varia- 
tions in the calculated physical properties is mainly due to 
the variations in the key structural parameters, rather than 
changes in the screening of the Coulomb interactions. The 
detail of the method is included in online material. 
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Method 

We use the combination of density functional theory 
and dynamical mean field theory (DFT+DMFT) [^l] as 
implemented in Ref. S2, which is based on the full- 
potential linear augmented plane wave method implemented 
in Wien2K [S3 , to carry out our first principles calcula- 
tions. The electronic charge is computed self-consistently 
on DFT+DMFT density matrix. The quantum impurity 
problem is solved by the continuous time quantum Monte 
Carlo method [94l using Slater form of the Coulomb 
repulsion in its fully rotational invariant form. 

We use the experimentally determined lattice structures, 
including the internal positions of the atoms, from Refs. [^S] 
- 9T9]. We use the paramagnetic tetragonal lattice struc- 
tures, and neglect the weak structural distortions. This dis- 
tortion has a very small effect on the size of the magnetic 
moment, proving that the magnetism has electronic rather 
than structural origin. 

We studied the paramagnetic phase of all compounds at 
the same temperature of T = 116 K, and magnetic states 
(SDW and DSDW) at T = 72.5 K. Our ab initio estimation 
for the Coulomb interaction U and Hund's coupling J in 
BaFesAss are U = 5.0 eV and J = 0.8 eV [S20 . We checked 
that the Hund's coupling, which is very weakly screened in 
solids, is very similar in other compounds, such as FeTe, 
where it increases for less than 5% compared to BaFe2As2. 
As shown in Ref. [ 32T\ , physical properties are not sensitive 
to small variation of Hubbard interaction [/, hence we fixed 
Coulomb interaction U and J to the same ab initio deter- 
mined values {U = 5.0 eV and J = 0.8 eV) across all studied 
compounds in the paramagnetic phase and DSDW phase, 
whereas a fixed U = 5.0 eV and a smaller fixed J = 0.7 eV 
are used in the SDW state in consistent with our previous 
calculations for BaFe2As2 in the SDW phase. [S21j 
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FIG. SI: Atomic histogram The atomic histogram of the 
Fe-Sd shell for (a) FeTe and (b) LaFeAsO in the paramagnetic 
state and magnetic states. The 1024 possible atomic configura- 
tions are sorted by the number of 3d electrons of the individual 
configuration. 



To analyze the character of the many body wave func- 
tion, it is instructive to project it to momentum and real 
space basis. The one electron spectral function contains the 
information of the overlap between the many body wave 
function and the plane wave basis. It is also instructive 
to project the wave function to the atomic basis on the 
iron site. This projection can be presented in the form 
of the atomic histogram. In any given period of time, an 
iron atom visits many states from the atomic basis on the 
timescale, which is proportional to the quasiparticle mass 
enhancement. The probability to find an iron atom in the 
solid in one of the atomic states, is called the atomic his- 



togram, and a typical example is provided in Fig. [ST] Here 
the atomic basis is constructed from the five 3d orbitals 
of an iron atom, which together with spin, span a Hilbert 
space of size 2^^ = 1024. We sort these states first accord- 
ing to their occupancy N = 0, 1, ...10, and within the same 
occupancy, we sort them according to their atomic energy. 
Due to large Hund's coupling the first (last) few states at 
given N are the high (low) spin states. In Fig. SI we clearly 
see the spikes in probability for the high spin states (at the 
beginning of the constant N interval). Consequently, the 
low spin states (at the end of the constant N interval) lose 
substantial weight. In the absence of Hund's coupling, the 
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high and the low spin states would be equally probable. 
We plot in Fig. [Sl] two representative histograms, for FeTe 
and LaFeAsO compounds. The two histograms are qualita- 
tively similar, nevertheless the differentiation between the 
high-spin states and the low spin states in FeTe is more am- 
plified. In the magnetic states such as SDW and DSDW 
states, the high spin atomic states gain even more weight. 



as seen in Fig. IST 
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FIG. S2: Fe 3d DOS Atomic-like Fe 3d DOS for FeTe con- 
trasted with actual Fe 3d DOS of LaFeAsO and FeTe computed 
by DFT+DMFT. 

The valence histogram of a Hund's metal is fundamentally 
different from that of an oxide. While only a few atomic 
states have a significant probability in an oxide, Hund's met- 
als visit a large number of atomic states over time, resulting 
in a dramatic (40%) reduction of the magnetic moment due 
to valence fiuct nations. A monovalent histogram with only 
the atomic ground state would give iron magnetic moment 
of 4 /iB- 

Another interesting feature of Hund's metals is that very 
large number of atomic states has finite probability. For 
comparison, in transition metal oxides or in heavy fermion 
materials with similar mass enhancement as in iron pnic- 
tides and chalcogenides, the atomic histogram would con- 
tain only a small number of states with significant prob- 
ability [922]. Since the Hund's rule coupling J is equal 
to 0.8 eV, the energy spread of atomic states at constant 
N = 5 or N = 6 is very large, of the order of 6 — 7eV. 
Because there are many atomic states with finite probabil- 
ity that contribute to the one electron spectral function, 
and because those states are extended over a wide energy 
range, the spectral function does not have a very well de- 
fined atomic like excitations. To demonstrate this effect, 
an atomic spectral function of Fe 3d 



S2 



we plot in Fig. 
orbitals, obtained from the corresponding atomic Green's 
function defined by 



Kn|4|m)|2(P„ + P„ 



(1) 



3d orbitals, and are atomic probabilities displayed in 
Fig. [Sl] Clearly, the atomic spectral weight is distributed 
over a very large energy range. For comparison, a typical 
heavy fermion would have one sharp peak (a delta function) 
below the Fermi level, and another peak above the Fermi 
level, i.e., a lower and an upper Hubbard band. [922] 

In Fig. [S2];a) we also show the fuh DFT+DMFT spectral 
function of the iron atom in the solid for FeTe and LaFeAsO. 
One can notice that these spectral functions have a sharp 
quasiparticle peak close to the Fermi level. Due to larger 
mass enhancement in FeTe, the quasiparticle peak in this 
compound is substantially smaller than in LaFeAsO. The 
rest of the spectral weight does not have a well defined 
Hubbard like bands, not because the rest of the spectra 
would be coherent, but because of the unusual atomic his- 
tograms of the Hund's metals. A small feature around —2 to 
— 1 eV is however noticeably enhanced in FeTe compared to 
LaFeAsO. This peak was identified in Ref. S23 as an atomic- 
like excitation, which is found in atomic spectral function 
at —2.2 eV, and is related to the excitation from atomic 
ground state of d^ to atomic ground state of d^. 

In the manuscript, we showed that one important fac- 
tor in determining the size of the magnetic moment is the 
quasiparticle mass enhancement. Clearly the heavier quasi- 
particles with smaller quasiparticle effective width are more 
prone to ordering. It is interesting to inspect also the " quasi- 
particle height", i.e., the value of the one-electron spectral 
function at the Fermi level. In Stoner theory, this value 
plays a crucial role in determining the critical tempera- 
ture and the size of the ordered moment. In Fig. [SSj^a) 
we show the value of the density of states at the Fermi level 
in the paramagnetic state as obtained by both DFT and 
DFT+DMFT. Clearly, large density of states at the Fermi 
level is more compatible with the small moment rather than 



large moment (shown in Fig. S3'b)), which disfavors Stoner 



theory for explanation of the trends in magnetic states 
across iron pnictides and chalcogenides. 

We also show in Fig. [SSj^b) the magnetic moment in the 
SDW and DSDW phases calculated by DFT with both the 
local spin density approximation (LSDA[924]) and gener- 
alized gradient approximation (GGA[325]) exchange corre- 
lation functionals. We also repeat the paramagnetic fiuc- 
tuating moment and the experimental static ordered mo- 
ments from the manuscript for better comparison. It is clear 



from Fig. S3 'b) that the DFT calculated magnetic moments 



roughly follows the trend of the fiuctuating moment in the 
PM state, but is very different from the static ordered mo- 
ment, as already pointed out by Ref. £26. 



Optical properties 

Now we turn to the plasma frequencies in the paramag- 
netic state of iron pnictide and chalchogenide compounds, 
shown in Fig. 



S4 



where n, m run over all atomic states, and a runs over Fe 



We show separately the in-plane and c- 
axis values, as obtained by both the DFT+DMFT and DFT 
calculations. We also plot the experimentally determined 
in-plane values from Refs. [S27 for Nai-^FeAs, [ 328] for 
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FIG. S3: DOS and magnetic moment: (a) Total density 
of states at the Fermi level in the PM phase computed by DFT 
and DFT+DMFT. (b) The magnetic moment calculated by DFT 
with both LSDA and GGA exchange-correlation functionals in 
both the SDW phase and DSDW phase. The fluctuating moment 
in the PM phase calculated by DFT+DMFT and the experimen- 
tal magnetic moment in the magnetic states which are shown in 
Figl(a) in the manuscript and reproduced here for easier com- 
parison. 



BaFe2As2 and SrFesAss, [929] for LaFeAsO, and [930] for 
LaFePO. The DFT+DMFT calculated in-plane plasma fre- 
quencies agree well with existing optical measurements, but 
are significantly reduced from the DFT values, showing the 
important of correlation effect. The extracted plasma fre- 
quencies in the DFT+DMFT calculation for FeTe are most 
strongly reduced from DFT values, and bear bigger error 
bars due to the fact that the scattering rate in FeTe is so 
large that there is no well defined Drude peak in the optical 
conductivity. The c-axis Drude response is not very differ- 
ent from the DFT prediction, which can be rationalized by 
the fact that mass enhancement for the orbital is the 
smallest. 



In Fig. S4 we also plot the inverse of the interlayer dis- 
tance c* of the Fe planes, multiplied by a constant factor, 
i.e., 12/c*. There is clearly a correlation between the inverse 



DMFTab * Exp. ab 
12/c* LDAab 
DMFT c LDA c 




Compounds 

FIG. S4: Plasma frequency. The PM in-plane plasma 
frequency Uab and out-of-plane plasma frequency Uc for var- 
ious iron pnictides and iron chalcogenides calculated by both 
DFT+DMFT and DFT. The experimental PM in-plane plasma 
frequencies are taken from Refs. 3271- 330l 



of the interlayer distance and in-plane Drude peak strength, 
suggesting that crystal structure again plays the key role in 
optical response. 



Orbital blocking mechanism 

As explained in the manuscript, the orbital blocking 
mechanism dramatically increases the strength of correla- 
tions in some multiorbital systems with large Hund's cou- 
pling, when the atomic ground state has large spin S. 

To understand this mechanism in terms of mathe- 
matical equations, it is useful to translate the mul- 
tiorbital system to the language of a self consistent 
Kondo effect [^ST, describing the formation of compos- 
ite quasiparticles via a Kondo Hamiltonian having a form 
^k,k'a,/3 ^oi,(3dl,^d(3a'C^k'(3,a'Cka,a plus the band Hamilto- 
nian of free electrons and the impurity Hund's term. The 
Kondo interaction is given by /q,,/3 = (IdJ^^d/^D+c.c. with H 
an effective atomic Hamiltonian and <> an average of the 
most probable atomic configurations [S31 . In the limit of 
no Hund's coupling, the Kondo coupling / is independent 
of orbital indices and the characteristic coherence scale is 
given by Tk cx exp(— l/(/A^p)) where = 25* + 1 is the or- 
bital degeneracy. In the limit of large Hund's coupling, the 
impurity is represented by a maximal spin 5, and the 1^,13 is 
diagonal, since the Hund's rule coupling jams electrons into 
the same spin and different orbital one particle states in the 
low energy many body state. As a result, the Kondo energy 
is dramatically reduced to Tk cx ex.p{—N/{Ip)) [S32 . We 
clearly see that the Hund's rule coupling has a strong effect 
in reducing the coherence scale. 

The crystal field splittings regulate the orbital selectivity 
of the mass enhancement via the orbital selective blocking 
mechanism. The more diagonal the effective Kondo cou- 
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pling, the stronger the correlations. There are two types 
of high spin states in atomic 3d^ configurations, namely 
eg^t2g^ and eg'^t2g^. When the low energy atomic con- 
figuration is eg^t2g^ (e^^t2^^), the t2g (eg) orbitals remain 
blocked, while eg {t2g) orbitals can mix. From this consid- 
eration it is clear that the crystal field environment which 
puts eg orbitals below t2g makes the eg^t2g^ primary config- 
uration, and effectively blocks the t2g orbitals, causing sub- 
stantially larger effective mass for t2g orbitals. This crystal 
field sequence is realized in most of compounds we stud- 
ied. The only exception is LaFePO, in which eg states are 
slightly above t2g, leading to slightly stronger correlations 
in eg orbitals. 



Kinetic frustration 




Compounds 



FIG. S5: Hopping parameters The hopping parameters of xy 
and yz orbitals to some selected orbitals of their nearest neighbor 
and next nearest neighbor Fe atoms obtained from tight-binding 
downfolding the DFT bands of the PM phase. 

To gain some insights into the low energy physics of iron 
pnictide and chalchogenide compounds, we also constructed 
a low energy effective tight-binding model. This postpro- 
cessing step is here used to understand why a relatively 
small change in crystal structure can lead to enormous dif- 
ference in mass enhancement. When the model contains 
both the iron d electrons and pnictide/chalchogenide p elec- 
trons, the hopping parameters do not show any anomalies or 
specific trend which could explain the variation of masses. 
On the other hand, the effective model which contains iron 
d electrons only, gives a clear signature that the hopping of 
the electrons on the iron xy orbital is severely impeded on 
the more correlated end, from LiFeAs towards FeTe. These 
hoppings are displayed in Fig. [S5 



From geometrical considerations, it is clear that direct 
iron dxy — dxy overlap between two neighboring atoms is 
always negative. This is due to the sign alternation of the 
wave function on the lobes of a dxy orbital. On the other 
hand, the indirect hopping through pnictogen atom p or- 



bitals such as Px^ is positive. While the indirect hopping 
through pnictogen p orbitals is larger than direct hopping 
for xz and yz orbitals, this contribution is comparable for 
the xy orbital. When pnictogen height increases, the indi- 
rect hopping decreases, and almost exactly cancels the direct 
hopping, resulting in negligible effective low energy xy — xy 



hopping, hence tl{xy^xy) in Fig. S5 becomes vanishingly 
small in the FeTe end. It is important to note that the band- 
width of the xy orbital does not change very dramatically, 
because the off-diagonal hopping tl{xy,xz) remains large 
(0.22-0.24 eV) even in correlated compounds. Neverthe- 
less, the most important diagonal hopping tl{xy^xy) dra- 
matically decrease together with the next nearest neighbor 
t2{xy^xy) and off diagonal next nearest neighbor t2{xy^z'^) 
hopping. Due to this kinetic frustration mechanism the di- 
agonal hopping is small in correlated compounds such as 
FeTe. Due to orbital blocking mechanism, the mixing of 
the orbitals is blocked, and hence the correlations can in- 
crease dramatically in the xy orbital. 
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